McGill/96-22 

SOFT INTERACTION OF HEAVY FERMIONS* 



C.S. LAM 

Department of Physics, McGill University 
Montreal, Quebec, Canada H3A 2T8 
E-mail: Lam@physics.mcgill.ca 



Abstract 

We discuss a formula for the sum of tree diagrams of n-bosons 
scattering from a heavy fermion. This formula can be considered as a 
generalization of the eikonal formula to include non-abelian vertices. 
It will be used to demonstrate the consistency of multi-meson-fermion 
scattering in the large- N c limit. 



1 Introduction 

This work is carried out in collaboration with Keh-Fei Liu of the University 
of Kentucky 1 . The central piece is a tree-level formula for the sum of the n! 
diagrams describing the scattering of n bosons from a heavy fermion. This 
formula is kinematical and combinatorial, valid for any coupling and internal 
quantum numbers. It relies on the presence of a heavy fermion, just like the 
heavy-quark effective theory. It will be used to demonstrate the consistency 
of meson-baryon scattering in large-iV c QCD 2 for an arbitrary number of 
mesons, under certain criteria. We are in the process of extending it to study 
the consistency of large- N c QCD in loops, baryon-baryon scattering, and 
heavy quarks. 

This formula bears some superficial similarity to formulas for soft-pion 
productions 3 , in that multiple commutators appear in both. However, the 
present formula does not rely on chiral symmetry, nor current algebra. As 
stressed before, it is kinematical and combinatorial, it is generally valid under 
no particular dynamical assumptions. Its validity depends on the smallness 
of all energy scales compared to the heavy-fermion mass M, and not some 
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Figure 1: A tree diagram. Solid and Wavey lines are respectively fermions 
and bosons. 

dynamical quantity like f n . If we need an analogy with known formulas, then 
it is much more appropriate to think of it as a non-abelian generalization of 
the eikonal formula 4 . 



2 Sum of Feynman Diagrams 

Fig. 1 depicts one of the n\ tree diagrams for the process 

F'ip') - F{p) + B.ih) + B 2 (k 2 ) + ■■■ B n (k n ) , (1) 

in which an initial fermion F' emits n bosons E>i to become a final fermion 
F. The momenta of the particles are enclosed in parentheses. If a boson 
in (1) is moved from right to left, we will get an inelastic boson- fermion 
scattering process. Both possibilities are considered using the notation of 
(1). The initial and the final fermions are assumed to have a heavy mass M: 
pi = ( p ')2 = m 2 . 

We shall use the notation T a = [tit 2 t 3 • • • t n ] to denote a tree graph ob- 
tained from Fig. 1 by permuting the n boson lines. The index a runs from 1 
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to n\, and each ti takes on a different value between 1 and n. The numbers 
between the square brackets indicate the order of the boson lines from final 
to initial states. In this language Fig. 1 is denoted by Ti = [123 ■ ■ - n]. 

One encounters these diagrams in calculating the emission of soft photons 
from a (comparatively) heavy fermion. The result, summed over the n\ dia- 
grams, is described by the familiar 'eikonal formula' 4 , according to which the 
n photons may be treated as being emitted independently from the heavy 
source. Note that such factorization does not occur in each Feynman diagram 
so the sum is vastly simpler than the individual diagrams. Since we need to 
obtain a non-abelian generalization of the eikonal formula, it is useful first 
to recall how the abelian case is obtained. 

The ith propagator of Fig. 1 is given by 

Si = ^-Jt^ > H=P + i2 k i ■ ( 2 ) 

lvl i Hi lt j = \ 

For other permuted diagrams kj should be replaced by k tj , but in every 
case we have qo — P and q n = p' . Let u = p-ki/M be the energy of the 
ith boson in the rest frame of the final fermion. Our central assumption 
is uji <C M, i.e., that the interaction is soft compared to the mass of the 
fermion. Then the numerator of propagator (2) may be approximated by 
M + 7p provided the intermediate masses Mj are not too different from M: 
A Mi = (Mf - M 2 )/2M < M. In that case (2) can be approximated by 

Si = PA , (3) 

where 

A = n , (4) 

Wi = Yl)=\Uj is the total boson energy preceeding that propagator, and 
V = |(1 + 7p/M) is a projection operator. With this approximation, the 
vertex factor 7 M can be replaced by the velocity v M = p^/M of the fermion, 
because 

V^V = v^V . (5) 

For the abelian eikonal formula, a single fermion without resonance is 
considered, so AMj = 0. The off-shell (for line p') photon-emission amplitude 
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from Fig. 1 is then 

A* [T 1 ] = f[(2eve(k l )D l0 ) , (6) 

i=l 

where the subscript is there to remind us of the condition AM; = 0, and 
the asterisk denotes 'off-shell'. Note that Wi in A contains a sum of meson 
energies so there is no factorization here. 

With a certain amount of combinatorial algebra, it can be shown 4 that 
the total off-shell emission amplitude is 

At = £A*[T a ] = f[(- e -^) . (7) 

o=l i=l \ P' K i ) 

This eventual factorization of the boson variables then allows the conclusion 
that the n photons are emitted independently. 

For on-shell amplitudes, the condition M 2 — (p') 2 = is equivalent to the 
energy-conservation condition 

■y n n 

TrP-E** = E W i = ° ( 8 ) 

m i=l i=l 

in the approximation uji <C M. The on-shell amplitude Aq is equal to the 
off-shell amplitude Aq with the nth propagator removed, so it is 

A ~(M 2 -p' 2 )A* = (9) 

A more careful calcluation shows that A is actually proportional to the 
products of 8{uji) 4 . Formulas (7) and (9) will be referred to as the abelian 
eikonal formula. 

We shall now proceed to consider the non-abelian generalization of this 
formula. In that case, the on-shell amplitude from Fig. 1 is 

A[Tj] = ( P X\V 1 S 1 V 2 S 2 V 3 ■ • • K-iSn-iKlj/A') (10) 

where V« is the ith vertex attached to the boson Bi, Si is the propagator in 
(2), and A, A' represent helicities and other internal quantum numbers. We 
allow in this non-abelian version multi-channel configurations where both 
Vi and AMj should be regarded as matrices connecting fermions of various 
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helicities and internal quantum numbers. These matrices do not commute 
and it is in this sense that the amplitude is non-abelian. Within the same 
heavy- fermion approximation, the numerator M + jqi in (2) may be replaced 
by its on-shell value and written as J2x, M A i (9i)'"A i (?i)- 111 this way (10) can 
be written as 

A[7\] = V 1 D 1 V 2 D 2 V 3 ■ ■ • K-iA*-iK (11) 

provided we take 

Vi = (fcAilVilft-iAi-i) (12) 

and Di defined in (4). For the permuted diagrams, V* and AMj are to be 
replaced by V ti and AM (i . 

Eq. (11) will be our starting point for further analyses. As discussed 
above, it is valid to the leading order of uJi/M. However, if we interpret 
A[Ti] as the pole-part of the amplitude, then (11) is exact provided the 
bosons are massless and the boson momenta are all parallel to one another 
( 'forward scattering ') . This is so because under that condition • kj — so 
the expression in (3) and (4) is exact. So is eq. (8). Moreover, by taking 
the pole part of the amplitude, the numerator of (2) should be taken on shell 
(the residue of the poles) so the vertices can be rigorously replaced by their 
helicity matrix elements as in (12). 

Given that, it is more satisfying to regard (11) in this way, as an exact 
expression for the pole part of the forward massless scattering amplitude, so 
that an exact non-abelian formula for the sum can be obtained. At the end, 
we might want to put in boson masses, consider non-forward scattering, and 
the non-pole part of the tree amplitude. These lead to corrections of order 
Ui/M, so the formula would still be valid to the zeroth order of this small 
parameter. 

To proceed further, let us first continue to ignore the possibility of reso- 
nances and put all AMj = 0. Eq. (11) then becomes 

MTi] = V[Ti]Dq[Ti] , (13) 

where 

V[T 1 ] = V 1 V 2 V 3 ---V n . 1 V n , (14) 

and 

n— 1 n—1 i 

dati = n a„ = n ^— ■ us) 

t=l i=l L 
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Figure 2: Same as Fig. 1 but with spurious bosons (dashed lines) put in to 
take care of resonance-mass corrections. 

The sum over all diagrams A = Sa=iA)pa] can then be obtained by a 
combinatorial argument 1 to be 

A = in -£'v MC [nn]D [Tln], (16) 

6=1 

with the following notations. [Tin] = [t' bl t' b2 ■ ■ -t' b n] is the subset of trees 
T a with the boson B n fixed at the end, and the sum in (16) is over the 
(n — 1)! permutations of the remaining boson lines. V MC stands for 'multiple 
commutator' 

V MC [Tln] = [V bl [V b2 [V b3 ■ ■ ■ [V bn _,,V n ] ...]]]. (17) 

Line n is singled out in (16) but we could have singled out any other line 
instead. Note that the derivation of (16) is not trivial, as can be seen from 
the fact that each multiple commutator contains 2™ _1 terms so (16) contains 
a total of (n — l)!2 n ~ 1 terms, whereas (15) where this comes from contains a 
total of only n\ terms. 

Finally let us see how the formula changes when resonances are included 
to make A Mi ^ 0. In that case, instead of D i0 , we should use the full Di in 
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(4). By using the expansion 

oo 

A = Ao E (-AM.Ao) 4 , (18) 

we can cast this amplitude into sums of amplitudes without resonances by 
introducing spurious bosons. Specifically, for a given i = Y,7=i ?h the ex- 
panded amplitude can be represented by Fig. 2, which differs from Fig. 1 by 
having I spurious bosons (indicated by dotted lines) added. These dotted 
lines must appear between the first and the last wavey lines, but otherwise 
they are to be inserted in all possible ways. The dotted lines carry zero en- 
ergy and a vertex factor U ti = -AM tj , where U is the original (wavey line) 
boson immediately to the right of the dotted line concerned. All dotted lines 
between two consecutive wavey lines carry the same zero energy and possess 
the same vertex, so they are to be thought of as identical particles. In the 
presence of resonances, the sum over all diagrams are again given by (16), 
but now modified to describe Fig. 2 instead of Fig. 1, and a sum over all 
£ must also be taken. The permutation in (16) is similarly modified to in- 
clude all permutation of non-identical bosons in Fig. 2. Inspite of the many 
sums present, the important thing is that it remains to be given by multiple 
commutators of the V^'s and the U^s. 

3 Large-iV c QCD 

We can use (16) to demonstrate the consistency of inelastic meson-baryon 
scattering amplitude in large- N c QCD. The requirment of consistency for 
n = 2 and 3 led to very interesting physical predictions 2 . For an arbitrary 
n, to our knowledge the consistency has not been demonstrated and we shall 
use (16) to show it. 

Let me review briefly the problem encountered in the limit of large-iV c . 2 
The baryon mass M in that limit is of order N c , so from (2)-(4) the baryon 
propagator is of order 1. However, the meson-baryon Yukawa coupling con- 
stant grows like V^c? so a diagram like Fig. 1 will be of order N^ 2 . On the 
other hand, the complete n-meson baryon amplitude is known to go down 
with n like iV c 1_ra//2 . There is then a discrepency of n — 1 powers of N c between 
an individual diagram and the complete amplitude. How can one effect such 
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a huge cancellation between the individual diagrams in the sum? This is the 
consistency problem mentioned above. The problem of course gets worse as 
n increases. 

For n = 2 and 3, a multiple-commutator formula similar to (16) had been 
used to demonstrate the consistency 2 . It proceeds by showing that though 
individual vertices Vi are of order \f~N~ c so the product of n vertices behave 
like iV"/ 2 , the multiple commutators of these vertices in fact behaves like 
iV c 1 ~ n / 2 , as desired. The same proof works for an arbitrary n as long as a 
formula like (16) is available to convert the sum of diagrams into multiple 
commutators. 

To understand why multiple commutators lead to such an enormous can- 
cellation, we have to go back to study how a large- N c baryon is made of and 
how it couples. A baryon in this framework is a color-singlet object com- 
posed of N c quarks in the S state. Different baryons arise from different spin 
and flavor excitations, but the wave function must be symmetric in the spin 
and flavor variables of the N c quarks. Spin is described by the rotational 
group SU{2)j so the spin of a baryon is represented by the Young tableau 
in Fig. 3(a). Each box in the tableau represents a quark. Those columns 
with two boxes are spin singlets, and those columns with one box are spin 
doublets symmetrical under permutation of these columns, so the spin of the 
baryon is J if there are 2 J single-box columns. Fig. 3(b) is for the flavor 
group. The tableau is identical to the one in 3(a) to enforce the symmetry of 
the wave function. This time each box represents the flavor state of a quark. 
Under isotopic-spin symmetry the flavor group is SU(2)p. The total isospin 
of the baryon is / = J because there are 2 J = 21 single-box columns in 3(b). 
Putting these together we conclude that the allowed baryon resonances have 

T _ T _ 1 1 3 Nc 

i o 2 , x, 2 , 12' 

On the other hand, if we assume a flavor symmetry of SU(3)p, then the 
two-box columns in 3(b) give rise to a 3 representation while the single-box 
columns give rise to a 3 representation. If there are p single-box columns 
and q two-box columns in 3(b), so that N c = p + 2q, then the SU(3)p 
representation of the baryon is (p, q) (where (1,1) is the octet, (3,0) is the 
decuplet, (2,2) is the 27-plet, etc). For large N c , one or both of (p,q) must 
be large, and we are then talking about SU(3)p representaions too large 
to appear in Nature. Consequently it is useless to talk about large N c as 
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(a) 



(b) 



Figure 3: Young tableaux for (a) SU(2)j and (b) SU(2) F . 

an approximation to Nature if we insist on SU(3)f symmetry as well as 
the identical-particle symmetry of the baryonic wave function. In that sense 
large N c hates SU (3) p, which makes the sucessful broken-5'?7(3) p predictions 
from the large- N c analyses 2 even more intriguing. To consider large iV c as 
a useful tool we must go back to SU(2)p symmetry. We may still include 
strange and other quarks, but in no way they should be regarded as identical 
to the up and down quarks. New quarks of each species in the baryon are to 
be represented by a new tableau in Fig. 3(a), with 2K single-box columns 
if that is the number of new quarks of that species. They carry a total spin 
of K. There will be no corresponding tableau for these new quarks in 3(b) 
because they are strong-isospin singlets. Consequently the total spin carried 
by the up and down quarks is still /, so the allowed spin J of the baryon 
is obtained by adding I and K vectorially, i.e., it ranges from \I + K\ to 



We are now ready to describe how the baryon matrix elements for the 
vertices are computed. Let q\q be free-quark creation and annihilation op- 
erators. Let us abbrviate a one-body (free) quark operator as 



where V is either the (spin) Pauli matrices <r\ the flavor matrices t a , their 
products, or 1. For two flavors, t a = r a are the Pauli matrices, and for 
three flavors, they are the Gell-Mann matrices t a = X a . The use of SU(3)p 
quark operators and matrices does not imply that the baryon wavefunction 
is SU(3)p symmetric. 



I -K\. 



jTq = {T} , 



(17) 
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Any operator O can always be written as sums of products of operators 
of the form (17), the only constraint being that they must have the same 
quantum numbers as O. A product of m one-body operators will be referred 
to as an m-body operator, and it can be shown 2 that an explicit factor 
F m = A c 1_m must accompany the presence of an m-body operator. We need 
to find out the A c -dependence of the baryon matrix elements of {T} and 
its products. We are interested only in baryons of low J and / which exist 
in Nature, hence most of the columns (of order N c ) in Fig. 3 are two-box 
columns. They constitute singlet spins and singlet isospins and contribute 
zero to matrix elements of T = a 1 and r a . For these T, the baryon matrix 
element 



comes only from the single-box columns so it is of order 1. For T = 1, or o"V a , 
or if it involves matrices connecting the up/down quarks to other quarks, the 
two-box columns do contribute so the matrix element (r) is generally of order 
N c , the order of the number of columns. Similar arguments show that the 
matrix element ({T} m } of an m-body operator is at most of order N™, so 
with the factor F m in front, the matrix element of every term of O is at most 
of order N c (assuming that there are no explicit factors of N c accompanying 
the definition of O). 

Though ({r} m ) is generally of order N c , it might be as suppressed as 
jyi-m y: a ^ ^ e r' s are made up of a l, s and r a 's. When this occurs important 
phenomenological consequences often follow 2 . For example, if all the m > 1 
-body operators are suppressed, then the dominant contribution comes from 
the one-body operator, but that is just treating the quarks inside the baryon 
as if they were free, so this gives rise to a large- A c justification of the quark 
model. Another example has to do with the mass of a baryon. It can be 
computed from the matrix elements of the operators {1}, {J}-{J}/A C , etc., 
so it is of the form Mi = aN c + bJ(J + l)/N c . The mass split AMj is not 
only small compared to M as demanded by the approximations needed to 
reach (16), it is actually of order A"" 1 and not just order 1. Like the Skyrme 
model, or the strong-coupling model, the baryonic resonances are rotational 
levels very closely packed. 

It is now easy to see why (16) leads to the enormous cancellation between 
individual Feynman diagrams needed for the consistency. From free-quark 




(18) 
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commutation relation, the commutator of two 1-body operators is again a 
1-body operator, and the commutator of an m-body operator with an wl- 
body operator is an (m + ml — l)-body operator. As far as matrix elements 
are concerned, this means that each time we commute we lose one power of 
N c , so all together (n — 1) powers of N c are lost by the multiple commutator 
in (16), and this is precisely the amount needed to match the discrepency 
between individual and the sum of meson-baryon scattering diagrams. 

If baryon resonances are included, (16) can still be used with modifica- 
tions, as explained. The argument above valid for AMj = will still be valid 
if the additional vertices C/j = — AMj are of order 1. As we saw above, this 
is actually of order A 7 "" 1 so there is no problem when baryon resonances are 
present either. 

In this way we prove the consistency of the large- N c amplitudes if (16) is 
used, with or without baryonic resonances. Eq. (16) is exact for the pole part 
of the massless forward scattering amplitude. When boson masses are taken 
into account and the forward scattering restriction is relaxed, eq. (16) is only 
approximate. As such cancellations arising from the multiple commutators 
are still present, but (16) alone is not sufficient to demonstrate the complete 
cancellation of the (n — 1) powers of N c . We shall not discuss this more 
difficult problem here. 
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